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ABSTRACT 

Sco X-2 (GX 349+2) is a low-mass X-ray binary and Z source. We have analysed 
~156 ks of Rossi X-ray Timing Explorer data, obtained in 1998 January, on this 
source. We investigated the fast-time variability as a function of position on the Z 
track. During these observations, Sco X-2 traced out the most extensive Z track ever 
reported from this object, making this the most comprehensive study thus far. We 
found the broad peaked flaring branch noise that is typical of Sco X-2, with a centroid 
frequency in the range 3.3-5.8 Hz. We also discovered low frequency noise, and a 
new peaked noise feature, with centroid frequencies in the range 5.4-7.6 Hz and 11- 
54 Hz, respectively. We discuss the phenomenology of these features, their relationship 
with the power spectral components found in other low-mass X-ray binaries, and the 
implications for current models. In particular, the low frequency noise we observed 
was strongest at intermediate energies, in contrast to the low frequency noise seen 
in other Z sources. We also detected very low frequency noise, and have calculated 
complex cross spectra between intensity and hardness. We found that the very low 
frequency noise is not entirely due to motion along the Z track. 

Key words: accretion, accretion disks - stars: binaries: close - stars: individual: 
Sco X-2; GX 349-1-2 - stars: neutron - X-rays: stars 



1 INTRODUCTION 

Sco X-2 is one of the persistently bright neutron star low- 
mass X-ray binaries known as Z sources. Z sources are 
so-named because, over time, they trace out a 'Z' shaped 
track in an X-ray colour-colour (CD) and hardness- intensity 
(HID) diagram. The top of the 'Z' is called the horizon- 
tal branch (HB), the diagonal is the normal branch (NB), 
and the bott om i s the flaring branch (FB) (Hasinger & 
van der Klis 1989). Inferred mass-accretion rate increases 
from the HB, through the NB, to the FB. The mass- 
accretion rate is thought to reach the Eddington limit at 
the bottom of the NB, and is super-Eddington in the FB. 
The five other traditional Z sources are: Sco X-1, GX 17-1-2, 
Cyg X-2, GX 340+0 and GX 5-1. A related class of ob- 
jects are the atoll sources; they exhibit a curve d trac k in a 
colour-colour diagram (Hasinger & van der Klis 



1.1 Z and atoll source fast-time variability 

Through the study of the X-ray fast-time variability of 
low-mass X-ray binaries, we are able to investigate mass- 
accretion processes; in particular, the fast-time variability is 
a signature of the flow o f ma tter near the compact object 
(see review van der Klis 200C). The typical power spectral 
properties of Z and atoll sources are correl ated with position 
in the CD (Hasinger & van der Klis L989| ; also van der Klis 



1995aD 



1989) 
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On the HB, Z sources typically exhibit six variabil- 
ity componen ts tha t are phenomenologically related (e.g., 
Homan et al. [2002| and references therein): a band-limited 
noise component with cut-off frequency in the range 2- 
10 Hz, called low frequency noise (LFN); a quasi-periodic 
oscillation (QPO) in the range 15-60 Hz, known as the hor- 
izontal branch oscillation (HBO); a broad sub-HBO peak 
located between the LFN and HBO, and which has been 
identified as a sub-harmonic of the HBO; a second harmonic 
to the HBO; and a pair of QPOs in the range 200-1130 Hz, 
called kHz QPOs. 
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On the NB and lower part of the FB t here is a QPO 
in the range 6 -20 H z (e.g 
& van der Klis Eoool 



van de r Kli s 
Homan et al. |2002|). On the NB 



19954 Dieters 
it 

is known as a normal branch oscillation (NBO), and its 
frequency is roughly constant near 6 Hz. With movement 
through the NB/FB vertex the mean frequency and fuU- 
width-at-half-maximum (FWHM) of the QPO rise rapidly. 
It is then known as a flaring branch oscillation (FBO), and 
becomes indistinguishable from the underlying continuum 
at about 10-20 per cent of the way up that branch. 

Some of the features seen in the HB are also sometimes 
seen on the NB and FB. In Sco X-1, van der Klis et al. 



(1997) found that the NBO appeare d to ' peak-up' out of 



the LFN. In GX 17+2, Homan et al. (2002) found that the 



NBO and LFN existed simultaneously, and were separate 
features. They also observed a peaked noise feature in the 
FB, in addition to an FBO, which they identified as LFN. 

A power-law noise component, called very low frequency 
noise ( VLFN) , exists throughout the Z track, and is stronger 
at higher m ass- acc retion rates an d at h igher energies (e.g., 
van der Klis 1995a). van der Klis (1991) suggested that it is 



produced by the variations in intensity that are associated 
with movement along the Z track. Dieters & van der Klis 
( 200c ^constructed power spectra of rank number (see Sec- 
tion 2.1), in the frequency range 3xl0~^-0.01 Hz, and also 
of intensity, in the range 8xl0~''-0.2 Hz. They found that 
the power-law indices from the rank number power spectra 
were consistent with those from the intensity power specta, 
although we note that the uncertainties on the rank num- 
ber power-law indices were quite large. If VLFN is due to 
movement along the Z track, then the speed of that move- 
ment should be correla ted with the VLFN strength. Dieters 
& van der Klis (I2OOOI) found this to be the case in Sco X- 
1, Homan et al. (2002) found that the correlation exists in 
some parts of the Z track of GX 17+2, and Wijnands et al. 
( [l997| ) found that the VLFN strength in Cyg X-2 did not 
vary monotonically with speed. 

Finally, a noise component with a cut-off frequency in 
the range 30-100 Hz, called high frequ ency ri oise (HFN), has 



1995a) 



also been observed (e.g., van der Klis 

At their lowest inferred mass-accretion rates atoll 
sources exhibit: a band-limited noise component, modelled 
using a broken power-law with a break frequency in the 
range 0.1-32 Hz, called (rather confusingly) high frequency 
noise; low frequency QPOs or broad peaks in the range 1- 
67 Hz; and kHz QPOs. Except for the kHz QPOs, these 
features are very similar to those seen in bla ck ho le candi- 
dates (BHCs) (e.g., Wijnands & van der Klis 1999 and ref- 



erences therein). In atoll souces VLFN dominates at higher 
mass-accretion rates. 

A correlation was found in Z sources and atoll sources 
between the frequencies of the kHz QPOs, and the HBO 
or 1-67 Hz QPO; the same correlation was found between 
features in t he pow er spectra of BHC s (Psa ltis, Belloni 
van der Klis 



1999a; also Psaltis et al. 1999b). Wijnands 



van der Klis (1999) discovered another correlation in atoll 



sources and BHCs, between the break frequency of the band- 
limited noise (the HFN in atoll sources), and the centroid 
frequency of a QPO or bump in the range 0.2-67 Hz (the 
1-67 Hz QPO in atoll sources). The same correlation was 
found to exist in Z sources between the LFN and sub-HBO 
peak, and also between the sub-HBO and HBO (Wijnands 



fc van der Klis 199!: ; also Homan et al. 200j ). v an der Klis 
([19944b), Belloni, Psaltis & van der Klis (^002]), and Sun- 



yaev & Revnivtsev (2000) have also made comparisons be- 
tween the power spectra of BHCs and neutron star LMXBs. 

There is currently much effort being put into develop- 
ing models that can explain these correlations. The discov- 
ery that the power spectral components in black hole candi- 
dates follow the same correlations seen in Z sources and atoll 
sources gives support to those models which do not depend 
on the presence or oth erwise of a solid surface or magnetic 
field (e.g., van der Klis |200C| ). 



1.2 Scorpius X-2 (GX 349+2) 

Except for the fact that Sco X-2 has never been observed 
in the HB, its hardness-intensity diagram looks nearly iden- 
tical to that of Sco X-1 (e.g., Schulz, Hasinger & Triimper 



198!;). Therefore, one might expect the fast-time variability 



of Sco X-2 to be similar to that seen i n Sco X-1. 

Ponman, Cooke & Stella (|l988|) analysed EXOSAT 



data and found broad peaked noise in the FB, with a cen- 
troid frequency, FWHM and fractional rms amplitude of 
roughly 5 Hz, 1 Hz and 3 per cent, respectively. Kuulkers 
& van der Klis (1998) analysed a small RXTE dataset and 
found that the peaked noise became narrower and moved 
to lower frequencie s with movement from the NB into the 
FB. O'Neill et al. (2001a) analysed Gmga data, and found 



peaked noise with a centroid frequency and FWHM in the 
range 4-7 Hz and 6-12 Hz, respectively. They found that it 
was strongest at about 10 per cent of the way up FB, and 
was present until halfway up the branch. They referred to 
the feature as flaring branch noise (FBN), and concluded 
that: FBN is not the same as atoll source high frequency 
noise; it is difficult to explain it with the model for typical 
norm al /flaring branch oscillations (Fortner, Lamb & Miller 
198£); and it resembles the peaked noise seen i n the F B of 



Cyg X-2 at low overall intensities. O'Neill et al. (2001a) also 
found very low frequency noise that became stronger with 
movement up the FB, as seen in ot her Z sources. 

Zhang, Strohmayer & Swank (1998) observed Sco X- 
2 with RXTE in 1998 January. During those observations 
Sco X-2 traced out the most extensive Z track ever seen from 
that object; most importantly, it was found higher up in the 
NB than ever before. A pair of kilohertz quasi-periodic os- 
cillations were detected at 712 Hz and 978 Hz, respectively, 
in the uppermost part of the observed NB. Agrawal & Bhat- 
tacharyya (2001) analysed ~40 ks of RXTE data obtained 
in 1998 September and October and also found an extended 
NB. They detected peaked noise similar to that found in 
previous observations, though its fractional rms amplitude 
did not vary smoothly, and the feature was detected as high 
as 80 per cent of the way up the FB. They found that the 
properties of the peaked noise were similar throughout the 
Z track, and thus concluded that the peaked noise seen in 
the NB was the same as that seen in the FB. A QPO, with a 
centroid frequency of 3.8 Hz, was detected (90 per cent con- 
fldence level) in the FB during a dip in the intensity that 
occurred after a flare. During the September and October 
observations the behaviour of the VLFN was very differ- 
ent to that seen previously: it was, on average, stronger in 
the NB than in the FB, and did not exhibit an increase in 
strength with movement up the FB. 
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On the basis of phenomenology, the peaked noise ob- 
served in the NB and FB of Sco X-2 i s unh kely to be an 
N/FBO (e.g., Agrawal & Bhattacharyya 2001 and references 
therein). First, the peaked noise seen in the NB of Sco X-2 
is much broader than an NBO. Second, as Sco X-2 moves 
from the NB, through the NB/FB vertex, and then up the 
FB, the centroid frequency and FWHM of the peaked noise 
decrease, and then stay roughly constant. This is in contrast 
to an N/FBO, of which, with the same movement along the 
track, the centroid frequency and FWHM increase rapidly. 
Third, the peaked noise in Sco X-2 has been seen as far as 
80 per cent of the way up the FB, while FBOs are only 
present in the low est 10 -20 per cent of that branch. 

Zhang et al. (|l99^ ) did not analyse the low frequency 
fast-time variability in detail. We, therefore, decided to re- 
analyse those data. Our primary objectives were: to look 
for the LFN and HBO that are expected to accompany kHz 
QPOs; to better understand the phenomenology of FBN; to 
study the properties of the kHz QPOs as a function of both 
position on the Z track and photon energy; and, to further 
investigate VLFN. A prelim inary a ccount of this work has 
been given by O'Neill et al. (|2QQlb|). 




10* 1.5x10* 2x10* 

ntensity (counts/s) 

Figure 1. Hardness-intensity diagram of Standard 2 data (16 s 
averages). The hardness is the counting rate ratio between the 
8.7-19.7 keV and 6.2-8.7 keV bands, and the intensity is the 
counting rate in the range 2.0-19.7 keV. The soUd line shows the 
spline used to calculate Sz. A representative error bar is shown. 



2 OBSERVATIONS AND ANALYSIS 

The Rossi X-r ay Ti ming Explorer {RXTE; Bradt, Roth- 
schild & Swank 1993) proportional counter array (PC A) ob- 
served Sco X-2 between 1 and 29 January 1998. The ~156 ks 
of data were recorded simultaneously in six different modes. 
The 'Standard 2' mode data had a time resolution of 16 s, 
and Sco X-2 was visible in the energy range 2.0-40 keV. 
Three 'single bit' modes and one 'event mode' were used, 
with a time resolution of 122 /j,s. The effective energy bound- 
aries of the single bit and event mode data were 2.0-5.1-6.5- 
8.7 keV and 8.7-40 keV, respectively. We only used data in 
which all five proportional counter units were operating. 

2.1 Construction of hardness-intensity diagram 
and power spectra 

A hardness-intensity diagram was constructed from the 
Standard 2 data using 16 s averages, and is shown in Fig. |l| 
The counting rates were corrected for detector deadtime 
and background. In the HID, the hardness is defined as the 
counting rate ratio between the energy ranges 8.7-19.7 keV 
and 6.2-8.7 keV, and the intensity is the counting rate in the 
range 2.0-19.7 keV. Fig. |l| shows the most complete Z track 
ever observed from Sco X-2 (see also Zhang et al. 1998). 



The counting rate at the top of the FB is a factor of 2.6 
greater than the intensity at the NB/FB vertex. In compar- 
ison, the factor was: 2. 2 in Gi nga data (1.1-16.7 keV, 64 s 
averages; O'Neill et al. |2001a| ); ~2.1 in t he EX OSAT data 
(1-10 keV, 64 s averages; Ponman et al. 1988); and ~2 in 
the 1998 September and October RXTE data (2-16 keV, 
256 s averages; Agrawal & Bhattacharyya ^001 ). We note 
that using lower energy bands and longer averages reduces 
the observed maximum intensity of the FB. 

Light curves were extracted from the single bit and 
event mode data, with a time resolution of 244 /is. We 
summed the energy resolved data to produce light curves 
with various energy ranges. We calculated power spectra 



from 16 s intervals in each light curve; each interval corre- 
sponded exactly to one point in the HID. To investigate the 
0.0625-300 Hz frequency range we used a logarithmic fre- 
quency rebinning scheme, while we used linear rebinning to 
investigate the kilohertz region. 

A spline was fitted through the HID, and rank numbers 
were calculated for each point. Rank number is a one dimen- 
sional measure of p osition in the Z track; it was introduced 
by Ha singe r et al. ( |l99C| ), and is symbolised as S z (Hertz 
et al. 1992; see also Dieters & van der Klis 2000). It is a 



measure of the inferred mass-accretion rate. We defined the 
NB/FB vertex as Sz=2, and the top of the FB as S^^S. 

Average power spectra were then constructed on the 
basis of rank number, and were normalised to give fractional 
rms amplitude squared per Hz. The error given for each 
mean Sz is the standard deviation of the rank numbers used 
in the average. 



2.2 Fitting of power spectra 

To investigate the low frequency fast-time variability, we fit- 
ted the resultant power spectra in the range 0.0625-300 Hz 
using: a power-law to represent very low frequency noise, 
Lorentzians as required to represent peaked noise, and a con- 
stant to represent the deadtime-affected white-noise level. 
We investigated the kilohertz quasi-periodic oscillations by 
fitting the resultant power spectra in the range 500-1200 Hz, 
using two Lorentzians and a constant. 

The VLFN fractional rms amplitude was measured by 
extrapolating the best-fitting power-law and integrating over 
the range 0.01-1 Hz. The addition of a Lorentzian com- 
ponent to a model was deemed necessary only when the 
probability of exceeding the xt of the best-fitting model 
was <0.001, and an F-test for the inclusion of the extra 
Lorentzian was significant at or above the 3a level. We fixed 
the centroid frequency and/or the FWHM of a Lorentzian 
component either when they could not be constrained, or 
when a fractional rms amplitude upper limit was being mea- 
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sured. Uncertainties in the best-fitting parameters were mea- 
sured using Ax^ ~ 1, and upper limits were measured using 
Ax^ = 2.71 (95 per cent confidence level). 

We corrected the best-fitting fractional rms amplitudes 
for differential deadtime and channel cross-talk as appro- 
priate (see Appendix A). Finally, we corrected the frac- 
tional rms amplit udes f or the background counting rate 
(e.g., van der Klis 1995b). In the low frequency power spec- 
tra, the reduction of the rms due to the binning of data was 
negligible; the binning correction factor applied to the rms 
of the kHz QPOs was 1.1 (van der Klis |l989|). 



2.3 Investigation of very low frequency noise 

To further investigate the VLFN, we constructed a HID from 
the single bit and event mode data using 32 s averages. We 
defined the hardness as the counting rate ratio between the 
energy ranges 8.7-40 keV and 6.5-8.7 keV, and the inten- 
sity as the counting rate in the range 2.0-40 keV. We calcu- 
lated complex Fourier spectra from both the hardness light 
curve and the deadtime corrected intensity light curve, using 
1 s time resolution and 32 s intervals. Each interval, along 
with its associated Fourier spectrum, corresponded exactly 
to one point in the HID. We used 1 s time resolution so 
we could directly investigate the fast-time variability in the 
time-domain. The Fourier spectra had a frequency range of 
0.03125-0.5 Hz. 

For each 32 s interval we calculated complex cross spec- 
tra between the hardness and intensity. We multiplied the 
complex Fourier amplitudes of the intensity, by the complex 
conjugate of the complex Fou rier a mplitudes of the hard- 
ness (e.g.. 



1987 ). At each frequency in 



van der Klis et al. 
the resultant cross spectrum there is a corresponding cross 
vector C, with a real and imaginary component. The phase 
of the cross vector is the phase difference between the vari- 
ations in hardness and the variations in intensity; a nega- 
tive phase means that hardness lags intensity. A phase of 0° 
means hardness and intensity are positively correlated, and 
a phase of 180° means they are anti-correlated. 

We then constructed average complex cross spectra. At 
each frequency, we averaged together the real and imagi- 
nary components from all cross spectra whose corresponding 
point in the HID was within a specified range of hardness 
and intensity. The mean amplitude and phase at each fre- 
quency was calculated from the mean real and imaginary 
components. The uncertainties in the mean real and imagi- 
nary components are the observed standard deviation in the 
mean of those components, and we determined the uncer- 
tainties in amplitude and phase by using error propagation. 
The phases at the Nyquist frequency are always either 0° or 
180° so they have no error bars. 

If, at a particular frequency, the hardness and intensity 
are uncorrelated, then the phases of the individual cross 
vectors corresponding to that frequency will be uniformly 
distributed between — vr and +n. In this situation, the real 
and imaginary components of the corresponding vector G 
in the average cross spectrum will have Gaussian distribu- 
tions with a mean of zero. The real component of G has 
a variance, cr^, that is equal to the variance of the imagi- 



2 



2M 



(1) 



where a^c\ and fj,ic\ are the observed standard deviation 
and mean, respectively, of the amplitudes of the set of cross 
vectors that were averaged together, and M is the number 
of vectors comprising that average. 

Therefore, in the absence of a correlation between inten- 
sity and hardness, the squared amplitude \Gf of the average 
cross vector is distributed such that the probability of ob- 
serving a |Gp greater t han a certain threshold level T is 
given by (Dieters et al. 



tion) 



where 



2000 



Dieters, private communica- 



Prob(|G[' > T) = exp(-r/|Go| 



|Go| 



|Go| 



2 , 2 

M 



(2) 



(3) 
(4) 

(5) 



Given an observed squared amplitude T, the significance of 
the detection of a correlation, expressed as a percentage, is 
given by 



significance = 100 x (1 - exp(-r/|Go| )) 



(6) 



nary component 
expression 



(jj . These may be calculated by using the 



The cross vector at the Nyquist frequency has only a 
real component, so it must be treated separately. The value 
Gp/(TR is distributed as with 1 degree of freedom, where 

'^l = Kc|+P|C|)/A^- 

For a comparison with the analytical significance lev- 
els, we also used a Monte Carlo simulation to determine the 
significance of the observed correlations. We randomised the 
phase difference between the hardness and intensity at each 
frequency, in each individual complex cross spectra, while 
maintaining the observed amplitudes. We then used those 
'random-phase' cross spectra to construct an average cross 
spectrum. We repeated this procedure 10000 times and ob- 
tained, at each frequency, the distribution of mean ampli- 
tudes that would be present if our data were actually uncor- 
related. With this method, the significance of a correlation 
is equal to the percentage of random-phase amplitudes that 
are less than the observed amplitude. 

We deemed the hardness and intensity to be correlated, 
at a particular frequency, if the significance obtained from 
the analytical method was at least 90 per cent. The phases 
of those cross vectors with a significant correlation are pre- 
sented in Table ^, along with the significance we obtained 
from each method. The results obtained from the Monte 
Carlo method were, on average, within 0.3 per cent signif- 
icance of those found from the analytical method (see Ta- 
ble |). 

We repeated the entire procedure twice. First, using a 
time resolution of 8 s with 256 s intervals, corresponding to 
the frequency range 0.00390625-0.0625 Hz, which allowed 
us to investigate the frequency range in which the power- 
law component dominates. Second, using a time resolution 
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of 0.125 s with 4 s intervals, corresponding to the frequency 
range 0.25-4.00 Hz, which allowed us to investigate the in- 
fluence of low frequency noise and flaring branch noise on 
our cross spectra. 



3 FAST-TIME VARIABILITY 

In Fig. ^ we present representative 0.0625-300 Hz power 
spectra from the NB (Fig. |^ top; Sz range 1.53-1.85) and 
FB (Fig. ^bottom; Sz range 2.03-3.0), and the various com- 
ponents we identified. We detected three features: very low 
frequency noise; low frequency noise/flaring branch noise; 
and a high frequency peaked noise component. 

We subdivided our data further and constructed thir- 
teen average power spectra covering the entire Z track. The 
best-fitting values from the thirteen average power spectra, 
in the range 0.0625-300 Hz, are presented in Table [ll 

In Table ^ we also present the I'max of each Lorentzian 
component. This can be calculated for both a Lorentzian and 
cut-off power-law, and it facilitates a comparison between 
the two forms. The !^max is the frequency at which there is a 
maximum in a power-density-times-frequency plot, and it is 
a useful measure of the characteris tic fre quency of a power 
spectral feature (e.g., Belloni et al. 2002). 



The kilohertz quasi-periodic oscillations were only de- 
tected at Sz=1.66. In Table ^ we present the best-fitting 
values of the kHz QPOs in the NB. We investigated the 
energy dependence of the kHz QPOs by fitting energy re- 
solved power spectra corresponding to Sz = 1.66. We fixed 
the centroid frequencies and FWHM at the values found 
from the 2.0-40 keV power spectrum. In the energy resolved 
spectra, we only detected kHz QPOs in the 5.1-40 keV 
range, with fractional rms amplitudes of 2.3±0.2 per cent 
and 2.1±0.2 per cent for the lower and upper peaks, respec- 
tively. In the energy range 2.0-8.7 keV, the upper limits of 
the rms for the lower and upper peaks were 2.1 per cent 
and 1.8 per cent, respectively. Therefore, the kHz QPOs 
in Sco X-2 are h ard, like those seen in other objects (e.g. 
van der Klis 



2000 



and references therein). Zhang et al. (199J ) 
fitted the kHz QPOs using Gaussians; we also fitted our 
QPOs using Gaussians and found best-fitting values consis- 
tent with their results. The /degrees of freedom (dof) was 
0.97/51 for the Gaussian model, compared with 0.93/51 for 
the Lorentzian model. 

To investigate the energy dependence of the 0.0625- 
300 Hz power spectral features we used energy resolved 
power spectra corresponding to the average spectra shown 
in Fig. bl We fixed the centroid frequencies and FWHM at 
the values found from the 2.0-40 keV spectra. The results 
are shown in Table ^. 

In the next three subsections we describe in more detail 
the properties of the LFN/FBN peak, the high frequency 
peak, and the VLFN. 



3.1 Low frequency noise/flaring branch noise 
peak 

In Fig. ^ we present the centroid frequency, FWHM, i^max, 
and fractional rms amplitude of the LFN/FBN peak (see 
Fig. ^ as a function of Sz . 



Very Low Frequency Noise N B 
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Figure 2. Representative power spectra (2.0-40 keV) from the 
NB (Sz range 1.53-1.85) and FB (Sz range 2.03-3.0). The white- 
noise levels have been subtracted and the solid lines show the 
best-fitting models. 



There is an increase in centroid frequency with move- 



ment down the NB, while the FWHM, 



and rms are 



consistent with being constant. The mean i/max in the NB 
was 8.2±0.2 Hz. With movement into the FB the centroid 
frequency decreases abruptly, and the rms increases sharply. 
There is no significant change in the FWHM and t'max. Over- 
all as the source moves up the FB the centroid frequency 
decreases while the FWHM increases. The I'max of the peak 
decreases as the object moves from Sz=2.05 to 2.15, and it 
then increases with further movement up the FB. We note 
that i^max varies less than the centroid frequency. The rms 
is strongest at Sz=2.05, and decreases until it becomes un- 
detectable in the top half of the FB. 

In the NB the peak is strongest in 6.5-8.7 keV band, 
and in the FB it is strongest in the highest energy band (Ta- 
ble |^. The overall phenomenology suggests that this peak 
is LFN in the NB, and then evolves into, or is replaced by, 
FBN as the source moves into the FB. This interpretation 
is discussed further in Section 



3.2 High frequency peak 
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Table 1. Best-fitting values from average power spectra (2.0-40 keV). 



VLFN LFN / FBN High Frequency Peak X?/dof 

Index rms Freq FWHM Vmax rms Freq FWHM i^max rms 
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(Hz) 




(per cent) 


(Hz) 


(Hz) 


(Hz) 


(per cent) 










1.66±0 


04 


1 


55±0 


07 


1.40±0.05 


5 


4±0 


2 


10.3±1.5 


7 


5±0 


5 


2.5±0.2 


22±11 


70±11 


41±7 


3.0±0.5 


1 


37 


/ 


76 


1.79±0 


04 


1 


37±0 


06 


1.33±0.04 


7 


6±0 


2 


7.3±0.8 


8 


4±0 


2 


2.3±0.2 


11±5 


57±11 


30±5 


3.3±0.3 


1 


25 


/ 


76 


1.91±0 


04 


1 


24±0 


04 


1.32±0.03 


7 


6±0 


2 


8.0±1.3 


8 


6±0 


4 


2.0±0.2 


14±4 


32±6 


21±3 


2.3±0.3 


1 


09 


/ 


76 


2.00±0 


02 


1 


56±0 


04 


1.83±0.05 


5 


8±0 


1 


11.8±0.3 


8 


3±0 


1 


4.06±0.04 


14" 


32" 


21" 


<1.0 





98 


/ 


79 


2.05±0 


01 


1 


68±0 


06 


1.99±0.08 


5 


5±0 


1 


11.7±0.3 


8 


0±0 


1 


4.33±0.05 


14" 


32" 


21" 


<1.7 


1 


14 


/ 


79 


2.09±0 


01 


1 


76±0 


06 


2.20±0.09 


5 


2±0 


1 


11.6±0.3 


7 


8±0 


1 


4.13±0.05 


14" 


32" 


21" 


<1.4 


1 


07 


/ 


79 


2.15±0 


02 


1 


74±0 


05 


2.35±0.08 


5 


4±0 


1 


10.3±0.3 


7 


5±0 


1 


3.72±0.04 


14" 


32" 


21" 


<1.6 


1 


03 


/ 


79 


2.21±0 


01 


1 


68±0 


04 


2.41±0.08 


5 


5±0 


1 


11.6±0.4 


8 


0±0 


2 


3.30±0.05 


14" 


32" 


21" 


<1.7 


1 


52 


/ 


79 


2.26±0 


02 


1 


72±0 


04 


2.78±0.09 


4 


6±0 


3 


14.3±0.7 


8 


5±0 


3 


3.01±0.07 


14" 


32" 


21" 


<1.7 





82 


/ 


79 


2.32±0 


02 


1 


71±0 


04 


2.99±0.09 


4 


4±0 


6 


16±1 


9 


1±0 


5 


2.52±0.09 


14" 


32" 


21" 


<1.4 


1 


34 


/ 


79 


2.41±0 


03 


1 


73±0 


04 


3.23±0.12 


3 


3±1 


5 


16±3 


8 


7±1 


5 


1.9±0.2 


14" 


32" 


21" 


<1.3 





96 


/ 


79 


2.55±0 


05 


1 


64±0 


03 


3.09±0.10 




3.3" 




11±2 


6 


4±0 


9 


1.3±0.1 


14" 


32" 


21" 


<1.1 


1 


14 


/ 


80 


2.77±0 


09 


1 


60±0 


02 


5.45±0.12 




3.3" 




11" 




6.4" 




<0.8 


14" 


32" 


21" 


<0.8 





97 


/ 


82 



"Fixed parameter. 



Table 2. Best-fitting kilohcrtz QPO values from NB average power spectra (2.0-40 keV). 







Lower kHz QPO 




Upper kHz QPO 




Xl/dof 




Freq (Hz) 


FWHM (Hz) rms (per cent) 


Freq (Hz) 


FWHM (Hz) rms (per cent) 




1.66±0.04 


715±12 


98±50 2.0±0.4 


985±7 


62±31 


1.6±0.3 


0.93 / 51 


1.79±0.04 


715" 


98" <1.4 


985" 


62" 


<1.1 


0.99 / 57 


1.91±0.04 


715" 


98" <1.2 


985" 


62" 


<0.4 


1.22 / 57 



"Fixed parameter. 



3.2.1 High frequency peak in the normal branch 

In the NB, the high frequency peak (see Fig. ^ top) had a 
centroid frequency in the range 11-22 Hz, and it exhibited 
no significant change in its centroid frequency or fractional 
rms amplitude (see Table |l|) . There is marginal evidence of 
a decrease in the FWHM with increasing rank number. In 
energy resolved power spectra in the NB we detected the 
peak only in the 8.7-40 keV band (see Table ^. 

In Fig. ^we present the 2.0-40 keV average power spec- 
trum from the uppermost part of the NB, with 82=1.66. We 
used power-density-times-frequency on the vertical axis be- 
cause this shows the power spectral features more clearly. A 
two-Lorentzian fit to the power spectrum (Fig. |^ top; Ta- 
ble 0) gave xS/dof=l. 37/76. The probability of exceeding 
this Xv is 0.0178, so the two-Lorentzian fit is formally sat- 
isfactory. But -x^ = 1.2i7 is still rather high. Moreover, there 
is a clear excess of power in the range 45-60 Hz, which is 
where we would expect to find a horizontal branch oscilla- 
tion. Therefore, we thought it was worthwhile to explore the 
option of fitting the NB data using three Lorentzians, even 
though the addition of another Lorentzian to the model is 
not statistically justified. 

We used a three-Lorentzian model to fit the 83 = 1.66 
power spectrum and obtained x^/dof=l. 16/73. An F-test 
indicated a 3.1o" improvement to the fit. The best-fitting 
three-Lorentzian model is shown in Fig. ^(bottom). We also 
fitted the power spectrum corresponding to 8z=1.79 in the 
same fashion and obtained x^/dof=l. 02/73. An F-test indi- 
cated a 3.5cr improvement to the fit. The best-fitting values 



from both of the three-Lorentzian fits are presented in Ta- 
ble |. 

We also fitted the 8z = 1.91 power spectrum with a three- 
Lorentzian model and obtained x^/dof=l. 06/73. An F-test 
indicated this to be a 1.4cr improvement. It is not surprising 
that, in this power spectrum, the three-Lorentzian model 
was not a significant improvement over the two-Lorentzian 
model. There is no obvious excess of power in the range 
45-60 Hz, in contrast to the excess power seen at 8z = 1.66. 
In Table W we show again, for comparison with the values 
from the three-Lorentzian fits at Sz = 1.66 and 1.79, the best- 
fitting values from a two-Lorentzian fit to the Sz = 1.91 power 
spectrum. 

Through comparison with the other Z sources, in partic- 
ular 8co X-1 and GX 17+2 (e.g., Wijnands fc va n der Klis 
199£ and references therein; Homan et al. 2002), we may 



identify the three Lorentzians as: LFN, a possible sub-HBO, 
and a pos sible HBO. A detailed comparison is given in 8ec- 



tion 4.2.2 



If we assume we have detected three peaked noise fea- 
tures at Sz=1.66 and 1.79, and only two features at Sz=1.91, 
then, as 8co X-2 moves down the NB, the peak representing 
the possible HBO disappears, while the possible sub-HBO 
remains strong (see Table With movement down the NB, 
the z/max of the possible sub-HBO is consistent with remain- 
ing constant. 

Even though an F-test indicated that, at 82 = 1.66, 
the three-Lorentzian model was a 3.1o- improvement over 
the two-Lorentzian model, we stress again that the two- 
Lorentzian model was formally a satisfactory fit. There- 
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Table 3. Best-fitting values from energy resolved power spectra from tlie NB and FB. 



Energy Range 




VLFN 


Peaked Noise rms (per cent) 


x2/dof 


(keV) 


Index 


rms (per cent) 


T,FN/FBN 


rTi crn Hrpn Pfi^iV 

X±ltlll J- i ± "OjtV 




Normal Branch 












2.0-40 


1.44±0.04 


1.35±0.03 


2.48±0.17 


3.2±0.4 


1.30 / 76 


2.0-5.1 


1.10±0.12 


0.85±0.03 


2.11±0.07 


<3.0 


1.50 / 81 


5.1-6.5 


1.55±0.15 


1.41±0.08 


3.20±0.08 


<3.6 


0.82 / 81 


6.5-8.7 


1.68±0.12 


1.74±0.10 


3.d1±0.u7 


<4.1 


1.44 / 81 


8.7-40 


1.52±0.06 


2.24±0.06 


2.74±0.13 


4.4±0.3 


1.08 / 80 


Flaring Branch 












2.0-40 


1.72±0.02 


3.25±0.05 


2.98±0.04 




1.28 / 76 


2.0-8.7 


1.68±0.02 


2.93±0.04 


2.83±0.01 


<1.0 


1.03 / 81 


5.1-40 


1.73±0.02 


3.89±0.06 


3.45±0.06 


1.3±0.1 


1.31 / 80 


2.0-5.1 


1.66±0.02 


2.42±0.04 


2.40±0.03 


<1.1 


1.00 / 81 


5.1-6.5 


1.69±0.03 


3.22±0.06 


3.14±0.04 


<1.1 


1.26 / 81 


6.5-8.7 


1.70±0.02 


3.74±0.07 


3.49±0.03 


<2.0 


0.95 / 81 


8.7-40 


1.77±0.02 


4.71±0.09 


3.78±0.03 


<2.0 


1.20 / 81 



fore, even though we have investigated the use of a three- 
Lorentzian model, we cannot conclude we have detected an 
HBO. Since we have not detected an HBO, we can measure 
an upper limit by fixing the centroid frequency and FWHM 
at the values found in, for example, Sco X-1. In that ob- 
ject, the centroid frequency was ~45 Hz, and the F WHM 



was in the range ~8-40 Hz (van der Klis et al. 1997 ). The 
2.0-40 keV upper limits of an HBO at S2=1.66, with the 
FWHM fixed at 10, 20, 30 and 40 Hz, were 1.0, 1.6, 1.9 and 
2.1 per cent, respectively. 



3.2.2 High frequency peak in the flaring branch 

We initially fitted the average FB spectrum (2.0-40 keV; see 
Fig. 1^ bottom) with a power-law and single Lorentzian, rep- 
resenting very low frequency noise and flaring branch noise, 
respectively, which gave x?/dof=l. 67/79. The probability 
of exceeding that xi is 0.00017, indicating an unsatisfac- 
tory fit at the 99.98 per cent significance level. Therefore, 
we added a second Lorentzian to the model and obtained 
xS/dof=l. 28/76. An F-test indicated this to be a 4.1a im- 
provement to the fit. The best-fitting centroid frequency, 
FWHM and fractional rms amplitude of the additional peak 
were 54^4^ Hz, 96I34 Hz, and I.OIq'j per cent, respectively. 
The best-fitting two-Lorentzian model is shown as a solid 
line in Fig. ^ bottom. 

We also detected the high frequency peak in the 5.1- 
40 keV power spectrum from the FB, with an rms of 
1.3±0.1 per cent (see Table The nondetection of this 
feature in the 2.0-8.7 keV band, with an upper limit of 
1 per cent, shows that it is genuinely hard. 



3.3 Very low frequency noise 

Very low frequency noise was detected in all regions of the 
Z track. In Fig. |^ we present the power-law index and frac- 
tional rms amplitude of the VLFN as a function of Sz. The 
NB/FB vertex seems to be a critical position. The index be- 
comes less steep as the source moves down the NB, and then 
abruptly becomes steeper again at the vertex. The index 
remains roughly constant with movement up the FB until 



Sz~2.5, after which it decreases slightly. In the FB, the mean 
value of the index is 1.70±0.02. The rms is consistent with 
remaining constant in the NB, and then increases abruptly 
at the vertex. The VLFN then continues to strengthen with 
movement up the FB, and is strongest in the highest energy 
band. 



3. 3. 1 VLFN and motion through the Z track 

We constructed an average power spectrum from the rank 
interval 1.87-1.88, which corresponds to the vertical part 
of the NB in the HID. In this region of the NB there is 
no associated change in intensity with motion along the Z 
track. Therefore, if VLFN is associated purely with move- 
ment through the track, then it should be absent in this 
power spectrum. We found VLFN with a power-law index 
and rms of 1.1±0.2 and l.liO.l per cent, respectively. These 
are consistent with the best-fitting values of the VLFN found 
in other parts of the NB (see Table |l|). 

The detection of VLFN in the vertical part of the NB 
suggests that VLFN is not due solely to movement along 
the track. However, it is possible that within each stretch 
of data over which Fourier transforms were calculated, the 
source moves a substantial distance away from the vertical 
region of the NB, by moving up and down the NB, or even 
into the FB, and thereby producing associated changes in 
intensity. In Fig. ^ (top), we present the HID constructed 
from 32 s averages of single bit and event mode data. In 
Fig. 1^ (bottom) we have plotted, in a HID, all of the data 
points in the vertical part of the NB (hardness interval 1.03- 
1.08; see Fig. ^ top) using a time resolution of 1 s. If there 
was substantial movement up and down the NB within each 
32 s interval, then we would expect the 1 s points to exhibit 
an arc- like distribution, reflecting the shape of the NB. We 
observed a vertical distribution, suggesting that the changes 
in intensity comprising the VLFN are not due purely to 
movement along the Z track. However, we caution that the 
error bars are large. 

A more sensitive way of investigating the correlation 
between hardness and intensity is to use cross spectra (see 
Section With phase measurements, we can determine if 
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Table 4. Best-fitting values from the three-Lorentzian fits to the Sz=1.66 and 1.79 power spectra, and the two-Lorentzian fit to the 
82=1.91 power spectrum (2.0-40 keV). 







Low Frequency Noise 




Possible sub-HBO 




Possible HBO 




Xl/dof 




Freq 


FWHM 


I'max rms 


Freq 


FWHM 


'-'max 


rms 


Freq 


FWHM 


'-'max 


rms 






(Hz) 


(Hz) 


(Hz) (per cent) 


(Hz) 


(Hz) 


(Hz) 


(per cent) 


(Hz) 


(Hz) 


(Hz) 


(per cent) 




1.66±0.04 


5.1±0.3 


8±2 


6.5±0.7 2.2±0.4 


13±5 


29±10 


19±5 


2.3±0.7 


48±2 


26±8 


50±2 


1.5±0.3 


1.16 / 73 


1.79±0.04 


7.3±0.2 


7.1±0.8 


8.1±0.3 2.4±0.2 


16±2 


16±6 


18±2 


1.9±0.5 


46±10 


59±17 


55±10 


2.1±0.4 


1.02 / 73 


1.91±0.04 


7.6±0.2 


8. Oil. 3 


8.6±0.4 2.1±0.2 


14±4 


32±6 


21±3 


2.4±0.3 


46" 


59" 


55" 


<1.5 


1.09 / 76 



"Fixed parameter. 



the variability at a particular frequency is not entirely due to 
movement along the Z track. In the FB and the upper part 
of the NB, where the track exhibits an increase in hardness 
with increasing intensity (see Fig. movement along the 
track must produce a phase of 0°. In the lower part of the 
NB, where the branch exhibits a decrease in hardness with 
increasing intensity, movement along the track must produce 
a phase of 180° . 

In Table we present the phases from the 256 s, 32 s 
and 4 s intervals. We only show phases for the frequencies at 
which there was a detection of a correlation between hard- 
ness and intensity at the 90 per cent confidence level. The 
upper NB corresponds to data with hardness >1.06, and 
the lower NB to data in the hardness range 1.0-1.06. For 
the FB, we used data with intensity >9000 counts s~^. We 
also present the significance of the detection of the correla- 
tion at each frequency, as measured from both the analyt- 
ical and Monte Carlo method, and rounded to the nearest 
one-hundredth of a per cent. Note that >3cr detections were 
generally only obtained in the FB. 

Since we wish to investigate VLFN, we need to deter- 
mine the frequencies at which VLFN is the primary influ- 
ence on the phases. Above ~2 Hz the variations in inten- 
sity are dominated by LFN/FBN, which is stronger than 
the VLFN at those frequencies; at lower frequencies the ef- 
fects of LFN/FBN may still be present. At frequencies where 
LFN/FBN dominates, the phases are characterised by non- 
zero values (see Table |5|); however, we cannot, of course, 
assume that every observation of a non-zero phase is due to 
LFN/FBN. If we observe, at a particular frequency, a phase 
that is consistent with being zero, then we can conclude 
that the phases observed at lower frequencies must be due 
to VLFN. Any non-zero phases observed at higher frequen- 
cies may be due to either VLFN or LFN/FBN. In the upper 
NB and FB, all phases in the range ~0.004-0.34 Hz and 
~0. 004-1. 2 Hz, respectively, are consistent with being zero. 
In the lower NB: there is a non-zero phase at ~0.004 Hz; the 
phases from ~0.016 to ~0.09 Hz are consistent with being 
zero; and, there are non-zero phases at frequencies in the 
range 0.1875-3.75 Hz. We shall, therefore, use 0.1875 Hz as 
our very rough estimate of the lowest frequency at which the 
phases are primarily influenced by LFN/FBN, and thereby 
conservatively restrict our study of VLFN to frequencies less 
than 0.1875 Hz. 

In the FB and upper NB, the phases are consistent with 
being 0°, indicating that the variations in intensity may be 
due entirely to movement along the Z track. We combined 
the 256 s interval FB data into a single frequency bin, with 



a mean frequency of ~0.033 Hz, and measured a phase of 
— 3?1±1?4, which is also consistent with zero. 

In the lower NB, the phase at ~0.004 Hz is consis- 
tent with being 180°; therefore, at this frequency, the in- 
tensity and hardness are negatively correlated. In the fre- 
quency range ~0. 016-0. 094 Hz, the phases are consistent 
with being 0°; therefore, in this frequency range, the in- 
tensity and hardness are positively correlated. Since the Z 
track in the lower NB exhibits a negative correlation between 
hardness and intensity, the movement giving rise to the vari- 
ations in intensity in the range ~0. 016-0. 094 Hz must have 
a cross-track component. And, since we observed a negative 
correlation between intensity and hardness at ~0.004 Hz, 
the cross-track component must become more prominent at 
higher frequencies. 

It is important to note again the limitations of our mea- 
surements. Only about half of our detections are at the >3cr 
level. This is particularly a problem in the lower NB, where 
we have only ~15 ks of data. 



4 DISCUSSION 



Using Zhang et al.'s (1998) RXTE data, we have carried 
out the most comprehensive investigation of the X-ray fast- 
time variability of Sco X-2 thus far, as a function of position 
on the Z track. Our primary discoveries are: the properties 
of the low frequency noise/flaring branch noise peak, and 
the very low frequency noise, exhibit abrupt changes with 
movement from the NB into the FB; there is a peaked noise 
feature with a centroid frequency in the range 11-54 Hz, 
which, in the NB, may contain either one or two components; 
and, very low frequency noise is not entirely due to motion 
along the Z track. 



4.1 Low frequency noise and flaring branch noise 

The peak we observed in the NB and FB, with a centroid 
frequency in the range 3.3-7.6 Hz, is not an N/FBO: NBOs 
are much narrower than the peak we observed; the centroid 
frequency of an FBO increases rapidly with movement up 
the FB, in contrast to our peak; and FBOs are only present 
in the lowest 10-20 per cent of the FB, while the feature in 
Sco X-2 is present until about halfway up the branch (see 
also Section |l.2| ) . Given the presence of kHz QPOs in the 
upper NB, and the fact that the peak has a z^max in the 
range 7.5-9.1 Hz, which is consistent with the LFN found in 
other Z sources (e.g., GX 17 +2, Homan et al. |2002| ; Sco X-1, 
Dieters & van der Klis 2000 ) , we argue that we have detected 
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9 



Frequency Phase (°) (per cent significance: analytical/Monte Carlo) 

(Hz) Upper NB Lower NB FB 



0.00390625 


6±7 


(100/100) 


174±19 


(94.56/94.80) 


-4±2 


(100/100) 


0.0078125 


8±15 


(99.95/99.99) 


- 


- 


-3±3 


(100/100) 


0.01171875 


-9±16 


(99.67/99.68) 


- 


- 


-2±3 


(100/100) 


0.015625 


-26±14 


(99.76/99.84) 


-57±26 


(91.40/91.96) 


-4±3 


(100/100) 


0.01953125 


-33±18 


(98.83/98.78) 


- 


- 


-1±4 


(100/100) 


0.0234375 


27±16 


(99.10/99.27) 


- 


- 


0±5 


(100/100) 


0.02734375 


6±29 


(91.28/91.74) 


- 


- 


-3±6 


(100/100) 


0.03125" 


38±25 


(91.67/91.43) 


- 


- 


9±7 


(100/100) 


0.03125*' 


20±8 


(100/100) 


-48±21 


(97.43/97.53) 


-1±2 


(100/100) 


0.03515625 


- 


- 


- 


- 


-2±8 


(100/100) 


0.0390625 


5±22 


(97.36/97.38) 


- 


- 


-12±8 


(100/100) 


0.04296875 


28±21 


(96.37/96.30) 


-15±25 


(90.62/90.78) 


8±9 


(100/100) 


0.046875 


- 


- 


- 


- 


-13it8 


(100/100) 


0.05078125 


14±22 


(91.91/92.10) 


- 


- 


-4±8 


(100/100) 


0.0546875 


21±21 


(94.90/94.77) 


- 


- 


3±9 


(100/100) 


0.05859375 


- 


- 


- 


- 


3±8 


(100/100) 


0.0625" 


- 


- 


- 


- 





(99.94/99.98) 


0.0625'' 


12±13 


(99.99/100) 


- 


- 


-7±4 


(100/100) 


0.09375 


0±19 


(99.11/99.18) 


-8±20 


(97.70/97.78) 


-8±6 


(100/100) 


0.125 


-16±20 


(97.38/97.37) 


- 


- 


-12±6 


(100/100) 


0.15625 


-53±20 


(98.43/98.41) 


- 


- 


-13±10 


(100/100) 


0.1875 


- 


- 


176±24 


(94.12/93.74) 


-13±12 


(100/100) 


0.21875 


-13±18 


(99.27/99.33) 






15±18 


(99.55/99.59) 


0.25'' 














0.25= 


-25±22 


(96.34/96.26) 


- 


- 


-5±12 


(100/100) 


0.28125 










-24±26 


(91.46/90.99) 


0.3125 


- 


- 


- 


- 


- 


- 


0.34375 


-15±23 


(96.38/96.06) 


112±22 


(94.62/94.50) 


- 


- 


0.375 


- 


- 


-177±28 


(90.55/90.73) 


6±21 


(97.43/97.82) 


0.4062 


- 


- 


- 


- 


28±20 


(98.03/97.82) 


0.4375 


- 


- 


- 


- 


- 


- 


0.46875 


- 


- 


- 


- 


- 




0.5'' 


- 


- 


- 


- 


- 


- 


O.S'^ 


- 


- 


- 


- 


- 


- 


0.75 
1.0 


102±16 


(99.78/99.76) 


- 


- 


- 


- 


1.25 
1.5 


; 


; 


; 




35±26 


(90.59/90.29) 


1.75 


-162±20 


(98.53/98.44) 










2.0 


-178±25 


(91.79/91.25) 










2.25 
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LFN in the NB. In the FB, the peak is similar to the FBN 
found earlier in Sco X-2 (e.g., Agrawal & Bhattacharyya 



2001 and references therein), and we identify it as such. 



LFN is a hard phenomenon (e.g., van der Klis 1995a) 



For example, Jonker et al. ( 2000 ) found that the fractional 
rms amplitude of the LFN in GX 340+0 was more than a fac- 
tor of three greater in the 13-60 keV band than it was in the 
2-5 keV band. In contrast, the LFN we observed is strongest 



in the 6.5-8.7 keV band. This difference suggests that maybe 
we have not observed LFN, since different production mech- 
anisms give rise to different energy de penden ces (e.g., com- 
pare HBOs with NBOs, van der Klis 1995a and references 



therein). However, the fractional rms amplitude spectrum 
of NBOs, for exa mple, is known to differ between objects 
(e.g., van der Klis 1995a and refe rence s therein; Dieters & 



van der Klis 20QC; Dieters et al. 200C) and LFN may be- 
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Figure 3. Low frequency noise/flaring branch noise centroid fre- 
quency, FWHM, fmaxi and fractional rms amplitude as a function 
of rank number (2.0-40 keV). The open circles indicate fixed val- 
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Figure 4. Power spectrum (2.0-40 keV) from the uppermost 
part of the NB (Sz=1.66), with Power Density X Frequency on the 
vertical axis. The best-fitting two- and three-Lorentzian models 
are shown in the top and bottom panels respectively. 



have similarly difFerent per source. This suggestion can be 
confirmed through a sensitive investigation of the Sz de- 
pendence of the LFN that we detected, and by making a 
comparison with other Z sources. 

In previous observations of Sco X-2 the properties of 
the peaked noise in the NB was found to be similar to that 
seen in the FB, suggesting that it was th e sam e feature in 
both branches (Agrawal & Bhattacharyya pOOl ). We argue, 
for two reasons, that FBN is difi'erent to LFN: the proper- 
ties of the LFN/FBN peak, as seen in Fig. ^, exhibit abrupt 
changes at the NB/FB vertex; and the energy dependence of 
the peak differs between NB and FB. We note that accord- 
ing to the 'unified model' (Psaltis, Lamb & Miller 1995), 
the inferred mass-accretion rates are expected to be sub- 
Eddington in the NB, and super-Eddington in the FB, sug- 
gesting that the mechanism producing the variability in the 
NB may be different to that in the FB. However, the prop- 
erties of the N/FBO, for example, also differ between NB 
and FB, with the chan ges o ccurring near the vertex (e.g., 
Dieters & van der Klis |200C| ), and the NBO and FBO are 
still thought to be related. Likewise, there remains the pos- 
sibility that LFN evolves into FBN. 

Our observations of FBN support the conclusions 
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Figure 5. VLFN index and fractional rms amplitude as a func- 
tion of rank number (2.0—40 keV). 



reached by O'Neill et al. 
Section 



1.2 



(2001c), which were detailed in 
In particular, we note again the similarity be- 
tween FBN and the broad peak seen in the FB of Cyg X-2 at 
low o verall intensities (Kuulkers, Wijnands & van der Klis 



1991: ; see also O'Neill et al. 2001a; Agrawal & Bhattacharyya 
01). The feature in Cyg X-2 was fitted with a cut-off 
power-law, with a J^max of 10±3 Hz. This is consistent with 
the FBN we observed. The peak in Cyg X-2 was not observed 
in the NB, with a fractional rms amplitude upper limit of 
1.0 per cent, compared with the FB rms of 2.7±0.2 per cent. 
An investigation of the Sz dependence of the Cyg X-2 peak 
will reveal whether or not it is the same phenomenon as 
FBN. Furthermore, a stringent upper limit on the rms of 
the peak in the NB would support our suggestion that FBN 
is not the same as LFN. 

We now compare the LFN and FBN we observed in 
Sco X-2, with the LFN obser ved in the NB and FB of 
GX 17-1-2 (Homan et al. |2002| ). The t^max of the LFN in 
the NB of Sco X-2 is generally a few hertz lower than that 
seen in the NB of GX 17-1-2, and it is consistent with the 
I'max observed in the HB. There is no abrupt decrease in the 
centroid frequency and increase in the fractional rms am- 
phtude of the LFN in GX 17-1-2, with movement from the 
NB into the FB, suggesting that it may not be the same as 
Sco X-2 LFN/FBN. However, Homan et al. (|oo|) stated 
that their best-fitting values of LFN in the lower NB and 




.5x10^ 
ntensity (counts/s) 

Figure 6. Hardness-intensity diagrams of single bit and event 
mode data. The hardness is the counting rate ratio between the 
8.7—40 keV and 6.5-8.7 keV bands, and the intensity is the count- 
ing rate in the range 2.0-40 keV. The top panel shows all data, 
using 32 s averages. The bottom panel shows NB data in the hard- 
ness range 1.03—1.08, indicated by the horizontal lines in the top 
panel, using 1 s averages. Representative error bars are shown. 



lower FB were probably influenced by the N/FBO that was 
present, so the true LFN phenomenology may have been 
masked. The overall behaviour of t'max in the FB is similar 
in both objects, and the decrease in the fractional rms am- 
plitude of the LFN in GX 17-1-2, with movement up the FB, 
is remarkably similar to FBN. In contrast to the FBN we ob- 
served, the LFN in GX 17+2 was detected as high as about 
85 per cent of the way up t he FB . But, we note again that 
Agrawal & Bhattacharyya (2001) reported the detection of 
peaked noise in the FB of Sco X-2 at about 80 per cent of 
the way up the FB. 

In addition to GX 17+2, a close examination of the LFN 
in Sco X-1 is also desirable. Sco X-1 has been observed by 
RXTE, with 2-60 keV counting rates of --10'^ counts s"\ 
and these data are current ly ava ilable in the archive (e.g., 
Yu, van der Klis & Jonker pooi| ). Of course, more data on 
Sco X-2 are also required before sufficiently detailed com- 
parisons can be made. 
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4.2 High frequency peak in the normal branch 

In the NB, we detected a high frequency peaked noise fea- 
ture, with a centroid frequency in the range 11-22 Hz. This 
feature contains power at frequencies where we would ex- 
pect to find high frequency noise, a horizontal branch oscil- 
latio n, a nd/or a sub- horizontal branch oscillation (see Sec- 
1.1). We showed that the NB power spectrum could 



We find that, acc ording to the correlation found by Wij- 
( 1999| ), and in comparison to GX 17-)- 2 



nands & van der Klis 



tion 



be fitted satisfactorily with a two-Lorentzian model, using 
one Lorentzian for the low frequ ency noise, and one for the 
high frequency peak (see Section 3.2). However, at Sz = 1.66, 



there is an excess of power in the frequency range where we 
would expect to find an HBO (see Fig. ^). The addition of 
a third Lorentzian component significantly improved the fit. 
Therefore, we shall consider both two- and three-Lorentzian 
models for fitting the NB data. 



J^.2.1 Two-Lorentzian model for normal branch data 

The high frequency peak, if considered as a single 
Lorentzian, resembles Z source HFN. The t'max we found is 
consistent with that seen in other Z sources (e.g., Hasinger 
& van der Klis 1989), and it was stronger at higher ener- 
gies, which is typical of HFN (e.g.. Dieters & van der Klis 
200C ). However, there is some doubt as to the origin of HFN: 



it may be a n inst rumental effect from EXOSAT (Berger & 
van der Klis Fggi). Homan et al. (|2002|) did not detect HFN 



in a large RXTE dataset on GX 17+2, in contrast to previ- 
ous observations made by EXOSAT. They suggested that in 
the HB in the EXOSAT observations, where the HFN was 
too strong to be an instrumental effect, the feature being fit- 
ted in power spectra may have actually been the sub-HBO 
peak. 

Is then the high frequency peak, taken as a single 
Lorentzian, a sub-HBO peak? At the top of the observed NB, 
the LFN we detected had a i^max of 7.5±0.5 Hz. The centroid 
frequency, FWHM and i^max of the high frequency peak were 
22±11 Hz, 70±11 Hz and 41±7 Hz, respectively. In compar- 
ison, when the LFN in GX 17-1-2 had a t-max of ~7.5 Hz, 
its sub-HBO peak had a centroid frequency, FWHM, and 
i^max of 25.2zfc l.3 H z, 21±7 Hz and 27±2 Hz, respectively 
(Homan et al. 2002). The high frequency peak in Sco X-2 is 
at the expected centroid frequency of a sub-HBO, but there 
are some marked differences: we would expect to clearly de- 
tect an HBO when a sub-HBO peak is present; the high 
frequency peak is much broader than a typical sub-HBO, 
which results in a much higher fmax; and, with movement 
down the NB, the fmax of the high frequency peak in Sco X-2 
decreased to 21±3, while in GX 17-1-2 the i^max of the sub- 
HBO peak increased to a maximum of ~50 Hz. Therefore, 
we conclude that the high frequency peak, if considered as 
a single Lorentzian, can not be identified as a typical mani- 
festation of a sub-HBO. 



4-2.2 Three-Lorentzian model for normal branch data 

We now consider the possibility that the high frequency peak 
consists of two Lorentzian features. At Sz = 1.66, the peak 
could be fitted using two Lorentzians, representing a possible 
sub-HBO and HBO. In the same part of the NB, the lower 
kHz QPO is at 715±12 Hz. 



(Homan et al. 2002), the frequencies of the LFN and possible 
sub-HBO peak we observed in Sco X-2 are consistent with 
expectations, assuming the presence of an HBO at 48 Hz. 
According to the correlation found by Psaltis et al. 



(1999a), when the lower kHz QPO is at 715 Hz the HBO is 



inferred to have a centroid frequency of about 60 Hz. This is 
somewhat higher than the 48 Hz centroid frequency of our 
possible HBO. An HBO centroid frequency of 48 Hz cor- 
responds to a l ower k Hz QPO freque ncy o f about 600 Hz 
(Psaltis et al. |l999ej ; Homan et al. |2002| ). The HBO is 
Sco X-1 has only been detected with centroid frequencies 
in the range ^ 40-50 Hz, and Sco X-1 also lies off the Psaltis 
et al. (1999a) correlation; for example, its HBO was found 
at ~45 Hz, even when the l ower kHz QPOs was at ~875 Hz 
(van der Klis et al. 1997). Therefore, we speculate that 
Sco X-2 is similar to Sco X-1 in this regard, and predict 
that future observations may reveal an HBO in the same 
frequency range as seen in Sco X-1. 

There remains a problem, however, with the interpre- 
tation of our i/niax~20 Hz feature as being a Z source sub- 
HBO. If we assume we have detected an HBO and sub- 
HBO in the NB power spectra at Sz=1.66 and 1.79, then, 
as Sco X-2 moves down the NB, the sub-HBO peak should 
disappear before the HBO does. This effect is clearly seen 
in GX 17-1-2, in which the sub-HBO peak becomes unde- 
tectable in the middle part of the NB while the H BO i s 
present all the way down to the FB (Homan et al. 2002). 
In contrast, we find that our possible HBO component dis- 
appears by Sz = 1.91, while the possible sub-HBO compo- 
nent remains strong. Furthermore, the i^max and centroid 
frequency of the possible sub-HBO at S2 = 1.79 and 1.91, 
are significantly lower than the maxim um v alues observed 
in the NB in GX 17-1-2 (Homan et al. |2002| ). However, we 
note that the uncertainties on our measurements do allow 
for an increase in the possible sub-HBO centroid frequency 
and i/max. Therefore, although problematic, we cannot rule 
out the possibility that we have observed a manifestation of 
a sub-HBO slightly different to that seen in GX 17-1-2. 

In summary, if three Lorentzians are required to fit the 
power spectrum from the upper NB, we conclude that: the 
frequency of the lower kHz QPO, relative to the centroid 
frequency of the HBO, is similar to that seen in Sco X-1 
(van der Klis et al. 1997); and the feature with j^max~20 Hz 
may be a sub-HBO peak, though it must be a different man- 
ifestation to that seen in GX 17-1-2. 



4.3 Comparison with models 

Since fiaring branch noise is present at super-Eddington in- 
ferred mass-accretion rates, where the accreting material 
and outgoing radiation must be separated (e.g.. Lamb 1991), 
an ad aptation of the photon bubble model (Klein et al. 
1996a) may prove useful in explaining it. The photon bub- 



ble model was originally applied to X-ray pulsars, in which 
locally super-Eddington accretion takes place onto the po- 
lar cap s of a strongly magnetised neutron star (Klein et al. 
1996a ). Photon bubbles are formed in the settling mound be- 
low the accretion shock. The bubbles coalesce to form larger 
bubbles, and then rise up and lose their photons through dif- 
fusion into the accretion shock. The result is quasi-periodic 
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variability in the X-ray intensity, the frequency of which is 
primarily dependent upon the time it takes for the bub- 
bles to coalesce. Modelling showed that QPOs could be pro- 
duced at frequencies from as low as 20 Hz to as high as 
12000 Hz, al ong wi th a simultan eous p ower-law component 
(Klein et al. [L996a| ; Klein et al. |l996b| ). It would be inter- 
esting to see whether or not a similar process, operating in 
either a radial inflow or in the accretion funnel onto the 
neutron star magnetic poles, could give rise to FBN. 

In the magnetospheric beat frequency mo del fo r hori- 
zontal branch oscillations (e.g., Alpar & Shaham 198!: ; Lamb 
et al. [1984 Shibazaki & Lamb [l987| ), the HBO is predicted 
to be accompanied by a noise component at lower frequen- 
cies, with both features expected to be hard. In the power 
spectra of Z sources it is possible to identify either the LFN 
or the sub-HBO pea k as t he noise component (e.g., Wij- 
nands & van der Klis 199!:). The LFN that we observed in 



the NB was strongest at intermediate energies, suggesting 
that, in the context of the magnetospheric beat frequency 
model, either LFN is not related to the HBO, or the fea- 
ture we observed is not classical LFN. We have argued in 
Section that we have observed Z source LFN. If we are 
correct, then the sub-HBO peak, rather than LFN, might 
be identified as the noise component that is expected to ac- 
company the HBO. 

In the relativistic precession model, the upper and lower 
kHz QPOs, and the HBO, are all produced by t he sam e blob 
orbiting in the accretion disc (Stella & Vietr i 1998; Stella 
& Vietri 



1999 



Stella, Vietri & Morsink [19991 ). The model 
makes clear predictions of the correlations between those 
frequencies. Unfortunately, since we did not definitely de- 
tect an HBO, and the kHz QPOs were only detected in one 
part of the Z track, we cannot test those predictions. And, 
our HBO upper limits are too high to draw any meaning- 
ful conclusion about its presence or otherwise, since weaker 
HBO s hav e been observed (e.g., in Sco X-1; van der Klis 
et al. 1997). If we assume that we have detected an HBO at 



48 Hz, with kHz QPOs at 715 Hz and 985 Hz, then the mass 
and spin frequency of the neutron star can be inferred from 
the model. Judging from fig. 1 in Stella & Vietri (1999), the 
inferred mass is ~2 Mq; the other Z sources have inferred 
masses in the range ~1.8-2.2 Mq (Stella & Vietri [1999| ). 
The mass of the neutron star in Cyg X-2 has been dynami- 
cally established as 1.78±0.23 Mq, wh ich is consistent with 
the inferred masses (Oro sz fc Kuulkers 1999). Judging from 
fig. 1(b) in Stella et al. (1999), the inferred spin frequency, 
for neutron star equation of state AU and an assumed mass 
of 1.95 Mq, is ~600 Hz; the other Z sources have an i nferre d 
spin frequency in the range ~600-900 Hz (Stella et al. 1999). 



In the so-called transition layer model (e.g., Titarchuk, 
Osherovich & Kuznetsov 1999| ), Z source LFN is the same 
phenomenon as atoll source HFN, and the sub-HBO peak is 
the same phenomenon as atoll source 1-67 Hz QPOs. Both 
of these phenomena are due to oscillations in the transition 
layer between the disc and the neutron star. The HBO, the 
HBO second harmonic, and the upper kHz QPO, are fre- 
quencies at which a blob oscillates, having been thrown out 
of the transition layer into the magnetosphere. The lower 
kHz QPO is the orb ital fr equency at the inner edge of the 
disc. Recently, Wu (2001) found good agreement between 
the predictions of the model and observations of Z sources 
and atoll sources. It is unclear how the fractional rms ampli- 



tude of LFN is predicted to vary with energy because the ac- 
tual mechanism that produces th at var iability has not been 
specified (Titarchuk & Osherovich 1999). Naturally any such 
mechanism that is proposed must account for the energy 
dependence we observed in Sco X-2. If the transition layer 
model is correct, we can measure the angle 5 between the 
plane of the accretion disc an d the equator of the magneto- 
sphere (e.g., Titarchuk et al. 1999). If we assume an HBO 



frequency of 48 Hz, and use our observed kHz QPO frequen- 
cies, then we infer 5=5?6±0?3, which is within the 5?4-6?4 



range found from the other Z sources (Wu 2001 



The sub-HBO seen in Z sources may be interpreted as 
the fundamental QPO frequency, with the HBO and its sec- 
ond harmonic being interpreted as the secon d and fourth 
harmo nics of the su b-HB O (Jonker et al. 2000). Jonker et al. 
(|00|) (also Jonker |200]] ) found a broad feature in GX 5-1 
at 1.5 times the HBO frequency, making it the third har- 
monic of the sub-HBO. They suggested that a warped disc 
with a two-fold symmetry can provide a mechanism which 
produces odd harmonics that are broader than the even har- 
monics, and that this mechanism could be used as part of 
the relativistic precession model. The FWHM of the possible 
sub-HBO and HBO we observed at Sz=1.66 were 29±10 Hz 
and 26±8 Hz, respectively. The uncertainties here are too 
large to test the warped disc model. 



4.4 Very \oyv frequency noise 

4-4- 1 The importance of studying VLFN 

If VLFN is due to changes in intensity that are associated 
with the movement of a source along its Z track, then a 
model for VLFN will provide at least part of the mechanism 
that gives rise to the Z track. 

For example, in the unified model (e.g., Psaltis et al. 



1995) the movement along the Z track is driven by changes 



in the mass-accretion rate which alter the physical condi- 
tions in the region of the neutron star and inner accretion 
disc. The mass-accretion rate changes from sub- to super- 
Eddington with movement from the NB into the FB. In this 
particular context, a model for VLFN will provide the mech- 
anism that produces those c hange s in the accretion rate. 

Recently, van der Klis (2001) suggested that position 
along the Z track may not be related to the actual mass- 
accretion rate. Rather, he proposed that rank number is 
related to the ratio rj between the instantaneous accretion 
rate through the disc, and the long-term average of that 
rate. In this model, movement to a higher Sz requires ei- 
ther: an increase in the accretion rate; or for the accretion 
rate to remain constant immediately following a decrease. 
Conversely, movement to a lower Sz occurs either: when the 
accretion rate decreases; or when it remains constant im- 
mediately following an increase. Again, if VLFN is due to 
movement along the track, then the details of this model 
must be such that rj exhibits a power-law variability. 

Furthermore, the abrupt changes we observed in the 
properties of the VLFN with movement through the NB /FB 
vertex suggests that any model for VLFN should also pro- 
vide information about the physical basis of that transition. 
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4.4.2 Observations of VLFN 

We observed an increase in the fractional rms amplitude of 
the VLFN with increasing Sz, and it was stronger at higher 
energies. The same behaviour w as seen in the Ginga obser- 
vations of Sco X-2 (O'NeiUet al. |2001aD , and also in Sco X-1 
(Dieters & van der Klis 2000). However, our results are very 



4.4-3 VLFN and movement through the Z track 



Dieters & van der Klis (2000) identified VLFN in Sco X-1 



different to those from the 1998 September and October 
RXTE data on Sco X-2, in which the VLFN behaved some- 
what erratically, and was, on average, stron ger in the NB 
than i n the FB (Agrawal & Bhattacharyya 2001 ). Homan 
et al. (^002j ) found that, in GX 17+2, the rms of the VLFN 
did not increase smoothly with movement up the FB; as 
the object moved away from the NB/FB vertex, the rms 
decreased at first, reaching a local minimum at about one- 
quarter of the way up the FB, and then increased smoothly 
with increasing Sz. However, the VLFN in GX 17-1-2 was 
still, on average, stronger in the FB than in the NB. The 
reason for the differences between our results and those from 
the 1998 Septembe r and October data is unclear. Agrawal 
& Bhattacharyya (2001) used 8 s intervals when calculat- 
ing their power spectra, which meant their lowest frequency 
was 0.125 Hz, compared to 0.0625 Hz in our spectra; this 
may have produced some differences between the two sets of 
results due to interference, during the fitting process, from 
the peaked noise. 

The power-law index of the VLFN we observed de- 
creased with movement down the NB, abruptly increased 
at the NB/FB vertex, and then remained roughly constant 
with movement up the FB, decreasing slightly in the up- 
per part of the FB; in the FB, the mean value of the index 
was 1.70±0.02. In contrast to this, the index of the VLFN 
found in the 1998 September and October RXTE data ex- 
hibited no dependence on Sz; in the FB the index had a 



mean value of 1.64±0.07 (Agrawal & Bhattacharyya 2001). 
In Sco X-1, the power-law index increased with movement 
down the NB, from 1.2 to 1.7, although, similar to our re- 
sults, was constant in the F B, with a mean of 1.76±0.04 
(Dieters & van der Klis |200C| ). In GX 17-1-2, the index grad- 



ually increased with movement down the NB and reached 
a maximum at about 10 per cent of the way up the FB; it 
then decreased grad ually with further movement up the FB 
(Homan et al. 2002). Of note, however, is that the power- 
law index in GX 17-1-2 reached a local maximum of ~1.6, at 
roughly three-quarters of the way down the NB, and it then 
subsequently decreased to reach a local minimum, of ~0.8, 



near the NB/FB vertex (see fig. 11 in Homan et al. 2002) 
A sim ilar effect has also been seen in GX 5—1 (Jonker et al 
2002), and we may have observed the same thing in Sco X-2 
albeit with a much broader local minimum. 



As discussed in the previous section, understanding 
VLFN is important in understanding the overall behaviour 
of Z sources. Therefore, the presence of vastly different 
VLFN phenomenologies in Sco X-2, during different obser- 
vations, is an important result that needs to be clarified. Wi- 
jnands et al. (1997) found that, in Cyg X-2, the properties 



of the VLFN varied as the overall intensity level changed; 
variations in other components of the fast-time variability, 
and in the X-ray spectral properties, also accompanied those 
changes. The differences in the VLFN phenomenology in 
Sco X-2 may be related to the same effect. 



with movement along the Z track, while Homan et al. (2002) 



found that, in GX 17+2, the speed of movement through 
the Z track was not correlated wit h VLF N strength in all 
parts of the track. Wijnands et al. (1997) similarly investi- 



gated the VLFN in Cyg X-2 at different overall intensities. 
They found that the speed of movement through the track 
was highest at the highest overall intensity level, while the 
VLFN was strongest at the intermediate level; they, there- 
fore, concluded that the VLFN in the NB could not solely 
be due to movement along the Z track. 

We detected VLFN in the vertical part of the NB. 
Therefore, if VLFN is produced solely by intensity changes 
that are associated with movement along the track, then, 
within each interval, there must significant movement along 
the track to regions above and below the vertical part of the 
NB, such that there are variations in intensity. The absence 
of an arc-like distribution of the 1 s time resolution points 
suggests there is no such movement. The 1 s time-domain 
data points have very large uncertainties in hardness, so we 
used cross spectral analysis to investigate the relationship 
between intensity and hardness. 

In the FB and upper NB the variations in the intensity 
and hardness are consistent with being in phase with each 
other. This is consistent with the variations in intensity be- 
ing due only to movement along the Z track, though it does 
not rule out a cross-track component. 

In the lower NB, we found that at ~0.004 Hz the vari- 
ations in intensity were negatively correlated with the vari- 
ations in hardness, which is consistent with those variations 
also being due to movement up and down that part of the 
NB. However, at higher frequencies, we found that the inten- 
sity was positively correlated with hardness, indicating that, 
at those frequencies, the variations in intensity must be asso- 
ciated with movement that is not directly along the Z track. 
And, the negative correlation at ~0.004 Hz suggests that 
that cross-track component of motion is more prominent at 
higher frequencies. We caution, however, that these conclu- 
sions must remain tentative due to the low significance, and 
large uncertainties, of our phase measurements. 



4.4.4 Models for VLFN 

A model has been proposed to explain power-law variabil- 
ity in blac k hol e candidates (e.g., Mineshige, Takeuchi 
Nishimora 



1994). Their model can produce power-law fiuc- 



tuations in the mass-accretion rate from the inner accrection 
disc, even when the rate of mass-injection into the inner part 
of the disc is random. The predicted X-ray power spectrum 
exhibits a power-law index of about 1.6, with a flattening 
at low frequencies. This model can possibly be applied to 
the unified model, thereby providing the mechanism to pro- 
duce changes in the accretion rate. Whether or not it can 
produce a large enough range in mass- accreti on rat e (0.7- 
1.02 times the Eddington limit; Psaltis et al. 1995| ) and a 
low enough break frequency (Sco X-1 is predicted to have a 
turno ver time-scale of several days; Dieters & van der Klis 
200c) needs to be tested. 

An alternative, yet similar, model for V LFN i s the 'drip- 
ping handrail' model (e.g.. Young & Scargle |l99q) , in which 
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the amount of matter present at the inner edge of the accre- 
tion disc varies both quasi- and aperiodically. If luminosity is 
considered to be proportional to the mass at the inner edge 
of the disc, then both a power-law variability and a QPO 
are produced. The frequency of the QPO is proportional to 
accretion rate, so it cannot be identifi ed as flaring branch 
noise (see also Dieters & van der Klis 2000). If the rate at 



which matter leaves the disc is considered, then the power- 
law component disappears (Steiman-Cameron et al. 1994) 



This model cannot, therefore, produce a power-law variabil- 
ity in mass-accretion rate, and so it cannot be applied to the 
unified model. 

Time dependent nuclear burning on the surface of the 
neutron star has also been proposed to produce VLFN, with 
the power- law index expecte d to increase with increasing 
mass- accretion rate (Bildsten 1995). Dieters & van der Klis 
(2000) noted that: the predicted variability has a character- 
istic time-scale, in contrast to the observed VLFN in Sco X- 
1, which is a power-law over time-scales ranging from 10 s 
to about 14 h; and the energy available from nuclear burn- 
ing is not sufficient to produce the observed fractional rms 
amplitude. Therefore, it is difficult for nuclear burning alone 
to account for VLFN. 

Our results provide further evidence that the variations 
in intensity comprising the very low frequency noise are asso- 
ciated with movement that is not directly along the Z track. 
A possible scenario is that changes in mass-accretion rate, or 
?7, produce motion in which there is a component both along, 
and across, the track, and that the cross-track component is 
more prominent at higher frequencies. Alternatively, changes 
in mass-accretion rate or r) might only produce movement 
along the track, with nuclear burning providing a cross-track 
component. Except for the transition from the NB into the 
FB, we observed a general decrease in power-law index with 
increasing inferred mass-accretion rate. In contrast, the nu- 
clear burning model predicts an increase in index with in- 
creasing mass-accretion rate. Therefore, if nuclear burning 
is found to be contributing significantly to VLFN, the 'r; 
model' for the Z-track is preferred, since it does not specify 
a particular mass-accretion rate for each pa rt of the track. 
As pointed out by Dieters & van der Klis (2QQC), the two 
types of variations would need to combine and produce the 
observed VLFN power-law. 



5 CONCLUSIONS 

We have carried out the most comprehensive study to date, 
of the X-ray fast-time variability of Sco X-2 as a function of 
position on the Z track. We found low frequency noise in the 
NB, and typical Sco X-2 flaring branch noise in the FB, with 
centroid frequencies in the range 3.3-5.8 Hz and 5.4-7.6 Hz, 
respectively. The LFN was strongest at intermediate photon 
energies, while FBN was hard. The various models seeking 
to provide an explanantion for LFN in Z sources must ac- 
count for the energy dependence we observed. With regard 
t o FBN , we confirm the conclusions reached by O'Neill et al. 
(2001a), and again note the similarity between FBN and the 



peaked noise seen th e FB of Cyg X-2 at low overall intensi- 
ties (Kuulkers et al. 199£). We suggest that an adaptation 
of the photon bubble model may po ssibly provide an expla- 
nation for FBN (Klein et al. |l996a|) . 



We found a new peaked noise feature in the NB and 
FB, with a centroid frequency in the range 11-54 Hz. This 
high frequency peak was stronger at higher energies. It could 
be fitted satisfactorily with a single Lorentzian, but, in the 
NB, it could also be fitted with two Lorentzian components, 
representing a possible HBO and sub-HBO peak. 

We found very low frequency noise throughout the Z 
track. It exhibited abrupt changes in the power-law index 
and fractional rms amplitude with movement from the NB 
into the FB. We calculated complex cross spectra between 
intensity and hardness, and discovered that VLFN is not 
entirely due to movement along the Z track. We speculate 
that VLFN may be due to: variations in mass-accretion rate 
or rj that produce motion in the HID with a component both 
parallel, and perpendicular, to the track; or a combination 
of nuclear burning on the neutron star, and variations in 
accretion or t). 
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APPENDIX A: CORRECTING RXTE DATA 
FOR DIFFERENTIAL DEADTIME AND 
CHANNEL CROSS-TALK 

In the case of RXTE, the observed counting rate per detec- 
tor /i, and incident counting rate A, can be divided into four 
componentsQ:'Good Xenon' counting rate /Xxe; the 'Remain- 
ing Count' the counting rate from the propane layer /Xp; 
and the 'Very Large Event' counting rate /Xv. 

The counting rate of Good Xenon events is the X-ray 
counting rate from the source and background. The Remain- 
ing Count is the counting rate of events that are detected 
simultaneously at more than one anode. These coincident 
events are primarily due to particles, and, in bright sources, 
may also due to 'genuine' photons. In each detector, there 
is a propane layer in front of the xenon layers which is used 
primarily for anti-coincidence, and may also be used to de - 
tect source photons in the range 1-3 keV (Bradt et al. 1993). 



The Very Large Event counting rate is from events that are 
above the upper discriminator; these events are mainly due 
to particles. 

Each component contributes to the detector deadtime. 
The deadtime per event from each component is: Txc ~ 
10 ^s, Tt ~ 10 fis, Tp ~ 10 fj,s, and Tv = 150 /is (tv is a 
commandable instrument parameter). 

The incident and observed counting rates are related by 



/i = /A (Al) 
where / is the deadtime correction factor, and is given by 

/ = 1 - TxcMxc — Tr/ir — Tp/ip - Tv^iv (A2) 

and 



(A3) 



1 + TxcAxc + '?"rAr + '^pAp + Tv Xv 

The actual fractional rms amplitude rxe of variations 
in AxG may be related to the observed fractional rms am- 
plitude Txe of variations in /ixc by the following expression 
(van der Khs |l989|) 



where 



C = 



C 



Axe df 



f dAx 



(A4) 



(A5) 



^ http:/ /rxte. gsfc.nasa.gov/docs/xte/recipes/pca_deadtime. html 
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2000 3000 4000 

Good Xenon Count (counts/s/deteotor) 

Figure Al. Deadtime corrected Remaining Count versus dead- 
time corrected Good Xenon counting rate. A representative er- 
ror bar is shown in the upper right. The solid line shows the 
best-fitting model fitted in the Good Xenon counting range 3000— 
4600 counts s~^ detector"^, and extrapolated to lower counting 
rates. 




2000 3000 
Good Xenon Count (counts/s/detector) 



4000 



Figure A2. Deadtime corrected propane layer counting rate ver- 
sus deadtime corrected Good Xenon counting rate. A representa- 
tive error bar is shown in the upper right. The solid line shows 
the best-fitting model. 



To determine C, we need to evaluate df/dAxc- There- 



fore, we also nee d to determine how A 
on AxG- In Figs. Al, A2, and 



Ap, and A 



A3 



respectively, versus A 
from Standard 1 mode usin g 16 s averages 

In Fig. 
Axe >3000 counts 



we present plots of Ar, 
The data presented are 



Al 



detector 



the region 
there is a roughly 
linear relationship between the Good Xenon count- 
ing rate and Remaining Count. In the region 
Axe <3000 counts s~^ detector"^, the lower envelope 
of the plot also follows the linear relationship observed 
above 3000 counts s~^ detector"^, suggesting that the rela- 
tionship holds for all Axe- The best-fitting linear coefficient, 
as fitted in the Axe range 3000-4600 counts s~^ detector"^, 
was 0.085. The large scatter of points lying above the linear 



2000 3000 4000 

Good Xenon Count (counts/s/deteotor) 

Figure A3. Deadtime corrected Very Large Event counting rate 
versus deadtime corrected Good Xenon counting rate. A repre- 
sentative error bar is shown in the upper right. 



relationship, in the region Axe <3000 counts s~^ detector"^, 
can be attributed to events from particles. This interpreta- 
tion is supported by the presence of a similar scattering of 
points in Fi g [A3| . 

In Fig. |A2|, there is a quadratic relationship between Ap 
and Axe- The best-fitting quadratic and linear co efficients 
were — 1.165e-5 and 0.175, respectively. In Fig. A3, there is 
no clear correlation between Av and Ax 



From the plots presented in Figs. Al 



A2 



and 



A3 



find 



dX 
dX 



0.085 



dA 
dAv 



2- = -(2.33e-5)Axe + 0.175 




rfAx 

Now we can evaluate df/dX^c and find C 

df _ -Txe [1 + 0.085 - (2.33e-5)Axe + 0.175] 



dAxc 
df 



dX 
And thus 

C : 



(1 + TxcAxc + Tr Ar + '^P'^p + TvAv) 

= -rxc[1.26 - (2.33e-5)Axe]./'' 



rxe[1.26 - (2.33e-5)Axc]Axe/| 



(A6) 
(A7) 
(A8) 

(A9) 
(AlO) 

(All) 



We used the above expression to correct for differen- 
tial deadtime and channel cross-talk in the standard manner 



(e.g., Lewin et al. 1992) 



This paper has been typeset from a T[5X/ M^HjX file prepared 
by the author. 
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